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Free motion around black holes with discs or rings: 
between integrability and chaos — I 
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ABSTRACT 

Geodesic dynamics is regular in the fields of isolated stationary black holes. However, due to 
the presence of unstable periodic orbits, it easily becomes chaotic under various perturbations. 
Here we examine what amount of stochasticity is induced in Schwarzschild space-time by a 
presence of an additional source. Following astrophysical motivation, we specifically consider 
thin rings or discs lying symmetrically around the hole, and describe the total field in terms 
of exact static and axially symmetric solutions of Einstein's equations. The growth of chaos 
in time-like geodesic motion is illustrated on Poincare sections, on time series of position or 
velocity and their Fourier spectra, and on time evolution of the orbital "latitudinal action". The 
results are discussed in dependence on the mass and position of the ring/disc and on geodesic 
parameters (energy and angular momentum). In the Introduction, we also add an overview of 
the Uterature. 
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1 INTRODUCTION 

It became evident no later than the 19th century that linear dynamical systems do not always represent adequately what happens in nature 
and society. It was also found that non-linear systems may behave in a complicated, "chaotic" way, even if governed by fully deterministic 
rules and effects. And since the non-linearity occurs in so simple cases as a set of three Newtonian gravity centres, linear systems in fact 
proved to be just marginal tips within a vast non-linear tangle. In the 20th century, mathematical background for the dynamical systems was 
notably supplied by the ergodic theory and by the study of differential equations (KAM theorem). However, it was only during the past half- 
century that the "theory of chaos" really established, focusing on simple non-linear dynamical systems whose orbits show highly sensitive 
dependence on initial conditions. Computers and numerical techniques allowed to follow and visualise long-term evolution of these systems 
and chaotic processes were soon detected in almost every discipline — the days of non-linear science have come (' Scottl2007h . 

It is for about 40 years now that the chaotic behaviour has been studied in general relativity. The amount of chaos can be expected higher 
in such an inherently non-linear theory, allowing the space-time dynamics itself to turn irregular. The latter has naturally been examined on 
cosmological soluti ons, mainly on Bianchi IX (Mixmaster) models. While the secular debate about chaoticity and the very na t ure of these 



universes continues j. . . : ColevlliooiJpav & Lehneill2004l:lBenini & Montanilliooj : ISoares & StuchilboOSi ; iHeinzle et dlbood ; 



Buzzi et al . 

2007; Andriopoulos & LeachI l2008l ; lHeinzle & Ugglj|20(m one of the later attractors is to assess the effect of the cosmological constant 



and/vs. that of a scalar field wi t hin the FLRW-model d ynamics Ijoras & Stuc"hill2003l : Iparaoni et allbood ; iHrvcvna & Szvdtowskilliooi : 
Lukes-Gerakopoulos et al.| [ 2008l : Macieiewski et al.ll2008l) . (We only give several more recent references here. The early ones can be found 



Hobill 



;rakop| 
et alj l 



1 1994) . for example.) 



The other type of dynamical system investigated in general relativity is the test motion in a given space-time. In this field, developed 
over the last 20 years, several obvious lines of research arise. First, it is interesting to start from some background where the (geodesic, or 
electro-geodesic) flow is regular and then consider various types of its "perturbations" in order to study whether and how much stochasticity is 
induced. In gravitational systems, most common ingredients are (i) how the background field and interaction are described (in a Newtonian, 
pseudo-Newtonian, post-Newtonian, post-Minkowskian or exact relativistic manner; and with or without "dark energy" contribution), (ii) 
how the test body is described (mainly whether it is point-like, without structure, or endowed with more multipoles), and (iii) whether 
and how the gravitational and EM radiation emissions and corresponding reactions on the body are taken into account. Such a study actually 
reveals whether the respective approximations (neglects) are admissible in judgements concerning the long-term evolution of a given physical 
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system. The "disentanglement" of the effects of various approximations may be problematic, however, because the latter may furthermore 
be interlaced intricately. 

To top it all, there is the "dynamics" of the employed numerical code. With how much (and what kind of) chaos does the program itself 
and round-off errors confuse the results? How are these "machine instabil ities" enlaced with those coming from true physics? Interestingly, 
all this need not lead to desperation according to the shadowing theorem dSauer et al ] |l997h which states that, for some systems (surely for 
hyperbolic ones), every numerically produced "pseudo-trajectory" is shadowed by a bunch of "tr ue trajectories" ( all of which star t from 
slightly different initial conditions in general) that remain close to it for a certain ("long") time (see Lai et al. i 19991) : Shi et al, 1 2001 ); Judd 
( 12008) for more recent contributions). Anyway, the above uncertainties remind that it is important to have more ways of studying the dynamics 
and more indicators of chaos with precisely known interrelations. They also point out the preference for invariant measures of chaos and the 
principal role of analytical methods (like the Melnikov method, often used in relativity). Let us mention, apropos this emphasis, an important 
result of iMottea ( 12003,) : contrary to the opinion prevailing until then, it was shown that the sign of Lyapunov exponents is coordinate 



indepe ndent under the conditions when the exponents can be meaningful indicators of chaos (see also the recent sequel by iMotter & Saal 
(2003)). 



One of natural backgrounds for the study of general relativistic test-particle dynamics are those corresponding to gravit ational waves, be- 
cause the latter have no Newtonian analogue. The waves may influence the particles in a chaotic manner, as already shown bv lVarvoglis & Papadopoulo^ 



19921) : they considered charged particles moving in a uniform magnetic field and affected by a linearly polarised weak gravitational w ave 



propa gating perpendicular to the magnetic field-lines. Similar conclusi on was obtained for geode sic motion in pp-waves CPodolskv & Veselvl 
1998h and recently extended to the waves described by Kundt metrics jPodolsk^ & Kofroij|2007b : chaos was detected there by analysing the 



dimension of boundaries between the sets of initial conditions ("basins") leading to different kinds of final outcomes. (Fractal boundaries 
evidence chaos.) 

However, it is doubly apropos to seek for chaos around the third principal "non-Newtonian" prediction of general relativity (besides 
dynamical cosmology and gravitational waves) — around black holes. Actually, in ancient Greek chaos meant a gaping bottomless void, 
where everything falls endlessly. . . 



2 DYNAMICS OF TEST MOTION IN BLACK-HOLE FIELDS 



Let us start a summary of previous results concerning black-hole fields by quotation from IComish & Frankell (19971): ". . . even the most 
pristine black-hole space-time harbours the seeds of chaos in the form of isolated unstable orbits. A small perturbation causes these unsta- 
ble orbits to break out and infect large regions of phase space. Note that the experience with Newtonian systems is very misleading. For 
example, the Kepler problem has more integrals of motion than are needed for integrability. Keplerian systems are thus impervious to small 
perturbations. In contrast, black hole space-times are at the edge of chaos, just waiting for the proverbial butterfly to flap its wings.'0 

The proverbial flap has been provided by various kinds of perturbations and their effect explored by several methods. The first pertur- 
bation we mention is on the side of the test particle: it consists in 



2.1 endowing the particle with spin, 

which deflects its motion from a geodesic and thus away from a sufficient number of conserved integrals. ISuzuki & Maedal i 1997 ) performed 
long-term integration of the Mathisson-Papapetrou equations (describing the evolution of spinning test particles in a pole-dipole approxi- 
mation) in a Schwarzschild field, in order to show, on Poincare sections, that the spin really leads to chaos. Later lSuzuki & Maeda 



I999I) 



analysed the imprints of particle dynamics on gravitational waves emitte d and learned that in a chaotic mode the emitted power increases, the 
frequency spectrum is "blurred" and the waveforms may also be changed. Kao & Cho ( 2005h confirmed the spin-induced chaos by calculating 
the Melnikov integral. This integral along the unperturbed homoclinic orbit from Poisson brackets of original and perturbation Hamiltonians 
reveals a "transversal distance" in phase space between stable and unstable manifolds emanating from a homoclinic orbit. Infinitely repeating 
zeros of the integral indicate that the stabl e and unstable manifolds are entangled, thus evidencing the occurrence of chaos. The same system 
was also studied bv lKovama et alj ( 120070 who were able to classify its chaotic motions into the "1/frequency" and "white-noise" types on 
the basis of the power spectrum of time series of the particle's z-position. The "1/frequency", resp. "white-noise" spectrum was detected for 
smaller, resp. larger values of the particle's spin, when the motion tends to be rather wea kly, resp. strongly chaotic. 

The dynamics of spinning particles was studied in the Kerr field, too. lHartll (l2003iJ!5) showed, by computing the Lyapunov exponents, 
that it may become chaotic, but only for the spin values that are too large to be astrophysically realistic, or even too l arge for the pole-dipole 
test-particle approximation to be valid. Rotation of the hole was noticed to rather enhance the instabilitv. iHanl 1 2008a ) found the effect of the 
centre's rotation just opposi te, by plotting Poincare s ections and mainly by computing a modified Lyapunov exponent — the "fast Lyapunov 



indicator". In the meantime. 



Kiuchi & Maedal J2004) checked whether chaos reflects on emitted gravitational waves. They concluded that the 



^ The role of uns t able c ir cular orbits and their re lation to homoclinic (separatrix) and zoom-whirl orbits were recently clarified, in the Kerr field, by 
iLevin & Perez- Gi3 j2009l) : |Perez-Giz & Levii] j2QQ9D . The stability of circular orbits (the null ones, in particular) was also rela ted to quasi-normal modes 
of black holes in the Schwarzschild(-de Sitter) and higher-dimensional (Schwarzschild-Tangheriini and Myers-Perry) space-times jCardoso et alj200^ . 
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spectram of waves from non-chaotic orbits contains only discrete characteristic frequencies, whereas that from chaotic orbits is continuous 
with finite widths (the waveforms did not seem to be altered significantly). 
Now let us turn to much wider range of options of 



2.2 perturbing the background. 



Like in the above case of the test body, one can "endow" the background (in a more or less arti ficial manner) with some e xtra multipoles. In 
the literature, their effects on dynamics was mostly inferred from Poincare surfaces of section. ^ Vieira & Leteliei 1996a ) thus showed how 
geodesies around a Schwarzschild black hole become chaotic under the perturbation by a quadrupole, a nd mainly octupole mom ents (in a 
Newtonian limit, the resulting field is an analog of the Henon-Heiles potential). Chaos was also detected i Letelier & Vieirairi997ah around a 
slowly rotating black hole superposed with a dipolar halo; particles counter-rotating with respect to the black hole turned out more unstable 
than the co-rotating ones. Relativistic , linear ised and Newtonian core-shell fields (Schwarzschild centre plus dipole, quadrupole and octupole) 
were compared by Vieira & Letelier 1 19991) . The dynamics proved to be sensitive to breaking of the reflection symmetry (disappearance 
of the equatorial plane) which enhanced/inhibited the chaos in oblate/prol ate fields. Anyway , the re lativistic dynamics was significantly 
more chaotic than its Newtonian counterpart. The latter was confirmed bv iGueron & Letelien ( 1200 ih for a black hole (vs. a point centre) 
superposed with a dipolar field. However, the system became even more chaotic (than the exact case) when the Schwarzschild-like centre 
was simulated by a pseudo-Newtonian (P aczyhski-Wiita) potential; this instability was amplified further by incorporating special relativistic 
equation of motion. IWu & Huang ( l2003h tested, on geodesies around the Schwarzschild hole with a dipolar shell, a new coordinate-free 
presc ription for computati on of Lyapunov exponents. A dipolar term was also employed as a perturbation of the R eissner-Nordstr0m space- 
time dChen & Wang|2003h ; the geodesic structure again turned chaotic. In the meantime. iGueron & Letelien ( 120021) superposed a quadrupole 
with a static or rotating black hole and observed that prolate deformations lead to chaos, whereas the oblate-system dynamics remained 
regular; in case of the rotat ing black hole, chaos prefere ntially occurred among counter-rotating orbits. Finally we mention the usage of 
the basin-boundary method: Ide Moura & Letelien ( l2000al) checked the escape of free test particles from a static black hole superposed with 
a dipole, quadrupole and octupole terms, and found null geodesies more regular than the time-like ones. This conclusion was reached by 
analysing the dimension of the boundaries of initial-condition "basins" that lead to different escape endpoints (in isolated gravitational 
systems, there are three possible outcomes usually: fall to the centre, escape to infinity and eternal orbiting). 

Rather than to "switch on" the multipol es individually, i t is also possible to incorporate them by turning to some more general exact 
stationary (usually also axisymmetric) fields. ISota et al. l Jl996l) devoted their study to several static axisymmetric space-times, namely those 
of the Zipoy-Vorhees class describing the fields of finite axial rods (including e.g. the Schwarzschild and Curzon metrics as special cases), 
and also to the system of A'^ Curzon-type point singularities distributed along the symmetry axis. They diagnosed geodesic chaos both 
on Poincare maps and Lyapunov exponents, but mainly tried to relate the known "reliable" means of recognising chaos (detection of the 
homoclinic tangle, in particular) with methods based on curvature properties of the space-time background. Specifically, they studied tidal- 
matrix eigenvalues and curvature of the fictitious space ob tained by energy-dependen t conformal mapping within the Newtonian approach, 
while the Weyl-tensor eigenvalues in the relativistic case. IVieira & Leteliei Jl996bh replied, however, by demonstrating that the criterion 
suggested in previous paper is neither necessary nor sufficient for the occurrence of chaos. They admitted there surely exist links between 
global dynamics and local (curvature) properties of a certain "configuration manifold" of the system, but concluded that "any local analysis, 
in effective or even physical spaces, is far from being sufficient to predict a global phenomenon like chaotic motion". (The "geometric" 
approaches have nevertheless made a conside rable progress since, see concluding remarks.) Two other, quite recent papers also illustrate 
differences that not so rarely occur in the field. iDubeibe et al J ( 120070 analysed the geodesic dynamics of the Tomimatsu-Sato space-time and 
its deformed generalisation provided by a particular vacuum case of the five-parametric solution of Manko et al. (considered e.g. as a possible 
description of the neutron-star field). They affirmed regular dynamics in the Tomimatsu-Sato (S = 2) solution, while in its generalisation 
they found c haos for the oblate case instead of the prolate one; this is in contrast to the results of lGueron & Letelien ( l2002h . Soon afterwards, 
Ha3j2008bl) re-examined the case of the Manko et al. space-time and came to the opposite conclusion (prolate field chaotic, whereas oblate 
field regular). 

Among the more general stationary axisymmetric black-hole backgrounds, the most important are those which might have some astro- 
physical relevance. These correspond to "the most symmetric" type of perturba tions — by additional axially and reflectionally symmetric 
bodies (namely a disc or a "halo"). Referring to the conditions in galactic centres, Vokrouhlickv & Karaj i 1 9981) computed Newtonian trajec- 
tories around a point centre surrounded by a thin disc and monitored how the long-term evolution of their param eters responds to the disc's 
gravitational field as well as to successive mechanical interactions switched on to simulate crossings of the disc. Saa & Venegerolesl ( ll999[) 
solved a similar problem, with an infinite homogeneous thin discH around a static bl ack hole (re lativistic version) or a point mass (Newtonian 
version). The relativistic system turned out to be more chaotic than the Newtonian. Saa i 2000l) then compared the results (namely Poincare 
sections) with the (Newtonian) case involving a thick disc (homogeneous in radial direction). Chaos were rather attenuated by the disc thick- 
ness. To the contrary, perturbations non- symmetric wi t h resp ect to the disc plane (dipole) turned out to strongly enhance the chaos. The same 
Newtonian system was later studied by iKiuchi et al.l ( 120071) on Lyapunov exponents and on time evolution of orbital parameters (and the 



^ Note that the potentials of such discs raise to infinity in perpendicular direction. 
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corresponding radial-motion power spectrum). The authors mainly showed, however, that the degree of chaoticity clearly reflects on gravita- 
tional waves g enerated (according to the quadrupole formula) by the orbiting particles, namely on their forms and energy spectra. We note 
in passing that Ramos-Caro et al. i 20081) analysed Newtonian motion in the fields of several members of the Kalnajs disc family (sole discs, 
without the central body) as a possible approximation of the motion of stars in spiral galaxies. They found chaos in the case of disc-crossing 
orbits (just induced by discontinuity of the field, no mechanical interaction is involved), whereas regular motion otherwise. IWu & Zhang 
1 20061) recognised, on Poincare sections and Lyapunov exponents as well, geodesic chaos in the relativistic field of a Schwarzschild black 
hole superposed with discs of th e "pecul iar" family constructed bv lLemos & Letelien ( 119941) from identical halves of two dipoles. Returning 
to a direct astrophysical setting: iKaras & S ubr (2007) proved the relevance of dynamics in the centre-disc fields for the long-term evolution 
of stellar clusters in galactic nuclei. The disc was found to bring the main "instability", whereas the spherical-cluster mean field and perihe- 
lion advance (included on post-Newtonian level) rather regularise the dynamics, the total effect still being the increase of probability of close 
encounters between stars and the centre (probably a supermassive black hole) by a factor of 10. 



Among the astrophysically important "perturbations", there is also counting in a magnetic field. In one of early Dapers. lKaras & Vokrouhlickv 

discussed the dynamics of charged particles in the Ernst space-time (Schwarzschild black hole immersed in a magnetic field) with 
the help of Poincare sections and Lyapunov exponents. Quite recently, Takahashi & Koyama ( 2009) plotted the Poincare sections through 
the dynamics of charged particles in the field of a Kerr black hole surrounded by the Petterson's test dipole magnetic field. They found that 
the dragging effects tend to suppress the instability [f] they explain this "regularisation" by the fact that for fast rotating holes the integral 
J / \/ geeu^u^ dr of the test particle's latitudinal four- velocity over one orbital cycle (of length I) becomes approximately c onserved. 

Another important perturbation is the one by gravitational waves. In another seminal paper. Isombelli & Calzetta i I992I) used the Mel- 
nikov method in order to detect the influenc e of a small period i c pertu rbation (corresponding to weak gravitational waves) on the geodesic 
structure of Schwarzschild space-time. Later Letelier & Vieiral ( 1997hl) explored, also by Melnikov method, the chaos in test motion in the 
Xanthopoulos solution for a Schwarzschild black hole perturbed by a gravitational wavey 

Another type of perturbation has already a long history in Newtonian celestial mechanics: perturbation of the test-particle motion in 
a given background by an additional distant body — the Hill's problem. .Moeckel ( 1992) tackled its relativistic version and showed, by 
calculating the Melnikov integral, that the third body p erturbs the two-body homo clinic orbit into a tangle of intersecting stable and unstable 
manifolds, thus implying that the system went chaotic. I Comish & Frankell 1 1997 ) analysed the geodesic structure of the equatorial plane of 
the extreme Reissner-Nordstr0m black hole, perturbed by an additional distant body. They detected chaos by examining the basin boundarie s 
and pointed out that it may e.g. lead to the increased production of gravitational waves. Chaos was also ascertained by|SetouetalJj2005), 
in a N ewtonian setting with the primary mass endowed with a potential simulating the Schwarzschild field. Similarly, ISteklain & Letelier 
i 2003) compared the Hill problem including the Paczyhski-Wiita pseudo-Newtonian (Schwarzschild-like) potential and compared it with the 
original Newtonian version. They applied Poincare maps, Lyapunov exponents and basin boundaries on specific systems of (Sun + Earth 
+ Moon) and (Milky Way -I- cluster -I- star) parameters and concluded that the pseudo-Newtonian systems are usually — but not always — 
more unstable than their Newtonian counterparts. Recently the authors extended the ana lysis to one of pseudo-Newtonian systems trying 
to incorporate dragging effects due to the rotation of the centre ( Steklain & Letelie3l2009h : Poincare maps, Lyapunov exponents and fractal 
escape method revealed that the orbits counter-rotating with respect to the centre are more unstable than the co-rotating ones. 

When the "third body" is heavy and nearby, it can no longer be treated as a perturbation and a two-centre (or multi-centre) problem 
arises. The non-rotating centres can stay in equilibrium if both (resp. all) bear charges equal to their masses and electrostatic repulsion just 
balances gravitational attraction. In relativity it means black holes of the extreme Reissner-Nordstr0m type; their multiple field is described by 
the Majumdar-Papapetrou class of exact solutions. Examination of the dynamics of free test particles in the field of two such fixed black holes 
by Co ntopoulos actually triggered the research in the whole area; a survey of the early papers was e.g. given bv lContopoulos & Papadaki 
1 19931). These showed that unlike the Newtonian dyna mics in the tw o-centre backgroun d, the re lativistic problem is strongly chaotic both for 
photons and massive particles. This was confirmed bv lDettmann et al. ( 1994) who also dDettmann et al,.1995.) discussed th e appropriateness 
of various tests of chaos (Lyapunov exponents, topological methods, Poincare sections) in general relativity on this problem. ! Yurtseveti J 1 995 ) 
showed that the motion of photons in the double-centre Majumdar-Papapetrou field can be translated as geodesic motion in Riemannian space 
described by the spatial part of the metric. This space has negative curvature, which is known to be a source of chaos. A similar geometric 
approach was tested (and compa red with information from Lyapunov exponents) on static multi-black-hole space-times with c ylindr ical 
symmetry by Szvdtowski (1997). Chaotic scattering of photons by doubl e Majum dar-Papapetrou centre was studied by LevinI ( 19991) . A 
thorough discussion of photon escape from the system was provided bv lAlonso et al. (2008). They identified a set of unstable periodic 
orbits that form a fractal repeller, calculated its main characteristics (Lyapunov exponents, Hausdorff dimension and escape rate, among 
others) and compared them with the results of numerical simulations. They also considered the limit case when one of the black holes may 
be considered as a small perturbati on of the other (it is light and far away). A very thorough treatment of the subject was also presented 
by IContopoulos & Harsoulal 1 2005 ). It was mainly focused on (fractal) structure of the basins of attraction of the individual black holes 
and of the escape basin, and on its dependence on ratio between masses of the holes. Several papers generalised or modified the problem. 



A result that can be understood as related to this one was obtained bv lLetelier & Vieir3 Jl998h : in a Taub-NUT space-time, the chaos induced by dipolar 
perturbation is attenuated by the NUT pai'ameter (which acts as a "gravitomagnetic monopole", i.e. as a source of dragging). 

The inherent iffiness of the delicate problems of nonlinear evolution may be illustrated by putting together two of the questions raised above: (i) how chaos 
in the motion of particles reflects on gravitational waves emitted by them and (ii) how incident gravitational waves perturb the dynamics of particles. 
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Howard & W ilkerson analysed Newtonian dynamics in the field of two centres with a finite dipole and discussed the nature of its 

chaoticity in detail. A guirreg abiria (1997. ) showe d that in the case of four ex treme centres, the scattering of charged test particles is chaotic 
both in Newtonian and in relativistic description. Cornish & Gibbona ( 1997 ) also took into account a dilaton field and observed transition 
between regular and cha otic motion for p hotons and for extremally charged particles. (They added a useful discussion on fractal, topological 
and curvature methods.) Ide Moura & Letelier ( 2000b) considered a slightly different (and more realistic) system — of two uncharged black 
holes, but far away from each other. The approximated motion of photons in this field shows chaotic features: the boundary between capture, 
orbital and escape basins is fractal, becoming a self-similar Cantor s et in the limit of infinite distance between the centres, with its box- 
counting dimension falling off logarithmically. Hanan & RadJ 1 2007 ) recently generalised the problem of multiple extreme black holes to 
higher-dimensional space-times (with special attention to the 5-dimensional case) and found qualitative similarity of chaos in test dynamics 
there. 

Finally, let us quote the variant of the orbital dynamics where the own field of the orbiting body is not negligible. The problem thus 
ceases to count among those with given background and involves dynamics of the Einstein equations themselves. Compact-binary evolution 
is a prominent astrophysical example of such a challenge. It in fact embodies, in a true tangle, all the knottiness both on the general-relativity 
side and on the si de of non-linear dyn amics. The chaotic evolu tion of compact binaries was discussed, on second post-Newtonian level 
including spins, bv lLevinI 1 2000l 2003 ) and by Wu & Xie (2008); in the last paper, the "fast Lyapunov indicat or" was computed in order to 
deduce the ro les of system parameters. (A great deal of information is also given in the 3rd PN treatment by Grossman & Levin 1 2009h .) 
Cornish & Levin (2003 ) also compared it (in the case with zero spins) with /e.s?-particle dynamics in a Schwarzschild background. For the 
chaotic, post-Newtonian system the Lyapunov time was sometimes found comparable to the timescale for dissipation due to gravitational 
waves, which would indicate possible abating of chaos due to waves (were the Lyapunov time longer, chaos would be damped). However, 
the conclusiveness of any post-Newtonian approximation in the issue of chaos (of course) remains a question (Levin 2006). 

In the present paper we study the dynamics of time-like geodesies in exact space-times describing the fields of static and axially 
symmetric black holes surrounded by thin discs or rings. We try to show, on Poincare sections, on time series of position or velocity and 
their Fourier spectra, and on behaviour of the "latitudinal-action", how the original, completely integrable Schwarzschild geodesic dynamics 
grow chaotic when the external source is being "turned on", and we check how the evolution of phase-space properties depend on the values 
of geodesic constants. 



3 STATIC BLACK HOLE WITH A RING OR DISC 



We are interested in time-like geodesic dynamics in the field of a black hole surrounded by a ring or a thin disc. Let us restrict to the static 
and axially symmetric case and to discs without charges and currents and without radial pressure. For such sources, the space-time metric 
can be put into the Weyl fornjf] 



2i/j,2 I 2 — 2i/ J ,2 . 2A — 2i//j z , , 'A 

e (it + p e d(p + e (dp + dz 



(1) 



in (Weyl) cylindrical-type coordinates {t,p,z,(j]), with the two unknown functions i/, A only depending on p and z. The "gravitational 
potential" v satisfies the Laplace equation, so it superposes linearly, while A is found by quadrature 

^ = / p{[{Kpf -{'^,^f]<ip + '2'^.pK^<iz} , (2) 

J axis 

computed along any line within the vacuum region. (The requirement that the space is locally Euclidean on the axis implies A = on its part 
lying outside the horizon.) 

3.1 Schwarzschild black hole 



The Schwarzschild field (of mass M) is described by 



1 di + d2 

2 ^ di + d2 + 2M 



2 4E 



where 

di,2 
E = did2 



4M^ 1 , r(r~2M) 

= - In — ^ — — 

2 E 



^p2 + =p Af )2 = T- - Af =F M cos 6, 
[{p^ + z^ + Nff - Az^M^]^''' = (r - 



(3) 
(4) 



MY 



cos^ e. 



The second expressions are in Schwarzschild coordinates (r, 9) which are related to the Weyl coordinates by 



* Geometrised units ai'e used in which c = 
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G = 1. Cosmological constant is set to zero. Index-posed commas mean paitial differentiation. 
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P= \/r{r - 2M) sxaO, z = (r-M)cos9, 
or, in the opposite direction, 

r-M= i(d2+di), McosO = ^{d2- di). 

Superpositions of tlie Schwarzsciiild blaclc liole witli an additional static and axially symmetric source are often being represented in 
Schwarzschild coordinates. Outside of the thin source with no radial pressure, the complete metric reads 

J 2 /, 2A/\ 2e,,2 , ^'^ , 2 , 2 -20, 2A,^2 , • 2^ ,,2n 

as = — 11 \ e at -\ ar + r e [e aO +sm Oacp ), (5) 

r 

where 0{r, 9) is the potential of the external source and A(r, 6) = \ — Asci™ (Aschw is given by ((U). In the present paper, the Bach-Weyl 
ring and several thin discs described by exact solutions of Einstein equations will be chosen as sources surrounding the Schwarzschild black 
hole in a reflectionally symmetric manner (they will lie in the "equatorial plane" of the hole). 



3.2 Bach-Weyl ring 

The thin-ring solution of Bach & Wevl 1 1922 ) is given by complete Ist-kind Legendre elliptic integral ^^(A;) = /J^^^ 



Abw 



2MK{k) 



nl2 



[{p + b){-K'^ + 4k''^Kk + ik'^k'^k^) - 4pk^k'^{k''^ + 2) 



where A4 is the ring's mass and 6 its Weyl radius, K 



(6) 
(7) 



^^,k'^=<M^,k^ = l-k'^ = ^,^ndh,2^^ipTbr + zK 



3.3 Inverted counter-rotating Morgan-Morgan discs 



A physically reasonable family of annular thin-disc solutions was obtained by iLemos & Letelieg j 19940 by inversion (with respect to the 
rim) of the finit e Morgan-Morgan discs interpreted in terms of two non-interacting streams of dust counter-orbiting on circular geodesies 
I Semerakll2003h . The Newtonian densities of inverted discs are (m = 1,2,...) 



w^"'\p) = 



(2m)! 7r2p3 



1-4 

p2 



-1/2 



(8) 



where A4 is mass and b is the Weyl inner radius. The corresponding potential is quite lucid in oblate spheroidal coordinates {x, y) which are 
related to the Weyl coordinates by 



= + 1)(1 - y^), z = bxy 

22'"+i(m!)2M 



n=0 ^2n 
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(10) 



where (for n < m) 

(m) _ (-l)"(4n + l)(2n)!(m + n)! 



C. 



(n!)2(m-n)!(2m + 2n+l)! 



and P2n, Q2n denote Legendre polynomials and Legendre functions of the second kind. Note that in computations it is convenient to express 
the latter as 



Q2„(iX) = -iP2„(iX) arccotX 



■P,-l(iX)P2n-,(iX). 



3.4 Annular discs with power-law density profile 

Another family of annular thin discs is the one with power-law radial profile of density i Semerakl|2004 ) 
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Figure 1. Radial dependence of the equatorial potential of a Schwarzschild black hole surrounded by a Bach-Weyl ring (6) on Schwarzschild radius r = lOAf 
(left), resp. by a power-law disc U2I13) (of m = 4, n = 10) with Schwarzschild inner radius r = lOM (right). The relative masses of the ring are M / M = 0, 
0.2, 0.4, . . . , 2.0 (the singular "valley" of course becomes more profound along this sequence), while the relative masses of the disc are M/M = 0, 0.5, 1.0, 
. . . , 5.0. It is seen that the disc is only scarcely able to develop its own potential minimum near the black hole, even though it is fairly concentrated (towards 
the lim) and its mass is chosen extremely large. Radius is in the units of M, potential is dimensionless. 



where m and n are natural numbers and {a)m = r(a + m)/T{a) is the Pochhammer symbol. The corresponding fields behave somewhat 
more regular at the inner rim than those of the inverted Morgan-Morgan discs. For even nlf| the potential can be expanded as 



j=0 



(2+22) b2, 
V " /m+l 



(12) 



in the "inner region" {^J~^^~^~^ < b), while as 



( rn , n ) 
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V " /m+l ^'^ ' 



E 



A:!(m-fc)!(p2 + z2)- 



(13) 



in the "outer region" (-\/ + > b). 



3.5 Complete fields 

The complete fields are given by sum of the potentials, u = uschv, + i'- The Schwarzschild singularity is the strongest possible source of 
gravitation, so it is not surprising that the additional sources can hardly affect its field considerably without being unrealistically heavy (and 
compact). Figure [T] illustrates this on the Bach-Weyl ring (one-dimensional source) and on one of the (rather concentrated) power-law discs 
(two-dimensional source). Extended (three-dimensional) sources would have even weaker immediate influence on the black-hole field. 



4 TIME-LIKE GEODESICS IN WEYL FIELDS 

Massive free test particles move along time-like geodesies. In Weyl fields these read, in the Weyl coordinates, 

du* A f E \ 2E , „ 

-r- = ^[ = — 5^ i/.ptt" + z^.zit" , (14) 

At At y ~gtt J e^" 

Au''' A { I \ 2e^''l 



At At \ g^^ J p- 
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\p{v_pu'' + v^^u'') - u''], (15) 



"^"^ - +^^{l-pv,,) + {v,,^\,)[{u'f~{u^f]+2{v,,-\,)u'u\ (16) 



At e^^ p^e 



(u,, - A,.) - {u'f] + 2{u,, - A,p)u''u^ (17) 



At e^^ p2g2A 

^ When n is odd, the potentials comprise logaiithmic terms. We will not consider such a case here, 
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Figure 2. Radial shape of the effective potential squared in the equatorial plane of a Schwarzschild black hole surrounded by a Bach-Weyl ring on 
Schwarzschild radius r = lOA/ (left column) and r = 30M (middle column) and by a power-law disc (of m = 4, n = 10) with Schwarzschild inner radius 
r = 15M (right column). Relative masses of the ring are Ai /M = 0, 0.1, 0.2, .... 1.0, while relative masses of the disc are M/M = 0, 0.2, 0.4, . . . , 2.0 
(the curves shift downwards with increasing mass of the additional source; the topmost curves correspond to a pure Schwarzschild black hole). From top to 
bottom row, the specific angular momentum £ is set to 0, 2\/3M, 4M and 5M, where M is the black-hole mass. Radius is in the units of M, the effective 
potential is dimensionless (per particle's rest mass). Note that display proportions are only kept fixed within the columns. Basic observations: (i) depending on 
parameters of the external source, several different configurations of circular orbits are possible; both types of these (unstable and stable) may lie either above 
or below the source; (ii) the ring being much "stronger" (singular) source than the disc, it always induces an unstable circular orbit below; (iii) there may exist 
two different unstable circular orbits (and a stable one in between) between the horizon and the ring/disc; a marginal case also exists when the stable and the 
"outer" unstable of these orbits coalesce into a marginally stable one. 
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Figure 3. Contour lines of the effective potential squared (TsJ in the meridional plane (r sin 9, r cos 9) of a Schwarzschild black hole surrounded by a Bach- 
Weyl ring. In the left column, the ring is on Schwarzschild radius r = 30Af and has mass M = 0, 0AM, O.SM (from top to bottom); the angular momentum 
is £ = 4M. In the right column, the ring is on Schwarzschild radius r = 2bM and has mass M = 2M, 20M; the angular momentum is £ = 0. (We have 
chosen such unrealistic values of the ring mass in order to show that even the singular ring has to be extremely heavy in order to attract locally as strongly as 
the hole.) The coordinates are in the units of the black-hole mass M. 



where gr^j^ is the metric tensor, the test-particle four- velocity, r is its proper time and 

£ = E/m = —ut = e^'^u , I = Ljra = u,f, = p^e~^''u'^ 

are constants of the motion (energy and azimuthal angular momentum per unit mass with respect to an observer staying at rest at spatial 
infinity). In the vacuum regions (which means everywhere when there are no extended sources), the derivatives of A can be expressed as 

from the Einstein field equations. (However, A itself must in general be found numerically at each point by integrating this gradient.) Both 
the metric functions are continuous everywhere outside the horizon, they only have a finite jump in normal derivatives on the disc which is 
given by the "Newtonian density" of the disc w(p) (the one standing in the Poisson equation), 

u,^{p, z 0^) = -i^,z{p, 2-^0") = 2ttw{p). 

Note that the equations for and it'^ do not contain A, yet the (t, (;/))-motion is not fully determined by u, being coupled to the meridional 
components u'', which depend on A explicitly. 

The additional-source effect on geodesic motion can be tentatively estimated on effective potential for radial and/or latitudinal motion. 
From the four-velocity normalisation and using the metric one has 



2\ 

e 



[{u'f + r{T - 2M){u''f] ^£^- (Voff)', where {V^sf = (^i _ ^) (^1 + ^fj^ ^ e'^ . (18) 



Above the black-hole horizon, the l.h. side is non-negative, so the particle energy £ has to be at least equal to VeS at a given place. 

It is favourable that the effective potential only depends on u (not on A). In the equatorial plane {z — 0), for example, the external 
gravitational potential reads 



2Mk('-^ 

Hp) = i.Bw(p) = , \ ^ (19) 



7r(p + fo) 

for the Bach-Weyl ring l|6ll, while JSemerakll2003l) 
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for the inverted m-th Morgan-Morgan disc iWl and ( Semerakll2004 ) 
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V " /m + l 

2 ;,2, 



below the disc rim 



[^'2,(0) 



\ " /m+l 



,2j 



above the disc rim 



(22) 



(23) 



for the "(n, m)-th" power-law-density disc (n is even) ( 112113b . Illustrations of the equatorial effective potential for the fields including Bach- 
Weyl rings and power-law discs are given in figure |2] It is seen that if the external source is sufficiently heavy and compact, there may exist 
more than one family of unstable and/or stable circular orbits — the orbits going around the black hole alone (staying below the ring/disc) and 
those skirting both sources. However, the equatorial motion being integrable, its potential provides little clue as to the chaoticity of general 
geodesies. 

The behaviour of the effective potential in the perpendicular, meridional plane can be sketched by drawing there its contour lines, 
(Veff 6; t) = const. In figure|3] this is done for the angular momenta (. = 4M and (. = Q. 
The potentials simplify most notably on the axis. That of the Bach-Weyl ring reduces just to 



M 



+ fe2 



while the inverted m-th Morgan-Morgan disc yields i Semeraklboosh 

22'"+i(m!)2x 



v{z) = 



E^2:'iQ2„ 



and the (n, m)-th power-law-density disc i Semeralj|2004 ) 

(1 + ^)™, M ^(-ly 



v{z) = 



+ fe2 Z-^ kn + 2 



2Fl(-, 



kn 



1 + 62/^2 



(24) 



(25) 



(26) 



where 2-F1 is the Gauss hypergeometric function. The pure axial effective potential (involving 1 = necessarily) is not very interesting, 
however, it always grows monotonously to zero when receding from the black hole. 

Let us remark that the requirement of stationarity (or even staticity) of self-gravitating sources typically leads to the occurrence of 
supporting stresses within them. As opposed to the Newton's theory, these stresses contribute to gravity in general relativity, which may 
result in peculiar structure of space-time near the sources. This is most noticeable in the vicinity of most singula r sources. Indeed, the 
direction-dependent eff ect of ring sources (even including apparent repulsive features) was noticed on geodesies by Semerak et al. 1999[) 
and further explored bv iD' Afonseca et alj ( 120050 . 



4.1 Astrophysical motions in black-hole & ring/disc fields 

Long-term ring/disc structures probably exist around black holes in X-ray binaries and in galactic nuclei. These accretion configurations are 
usually modelled within test approximation, i.e. their own gravitational effect is neglected. In a "relativistic" region close to the compact 
centre, this is almost always apropos as far as intensity is concerned, but the additional matter may alter higher derivatives of the field which 
— vice versa — determine stability of its configuration. Actually, the most important parameters of the accretion disc like its position and 
extension may thus be very sensitive to the precise shape of the field which it itself generates. (This issue nicely illustrates the frequently 
stressed feedback between matter and its field.) 

In stellar binaries, however, the accretion structures are supposed to be rather light and the whole accretion system is quite strongly (and 
non-stationarily) perturbed by the binary companion. Hence, the study of self-gravitating rings or discs — and of the respective perturbed 
dynamics — seems to have more relevance for the galactic nuclei. There the black holes are surrounded, in addition to inner accretion 
discs, by central star cluster and by massive cold tori on larger radius. The behaviour of central stars e.g. in our Galaxy is not yet fully 
understood and it is natural to ask about the influence of the other elements on their dynamics. In the first approximation, the torus can be 
described by a potential of a thin ring, the inner disc by that of an infinitely thin disc and, possibly, the cluster by an additional spherical 
term. In the present paper, we will thus be interested in motions of free test particles ("stars") below and above the radius of the ring/disc 
and monitor, in particular, whether their orbital parameters cannot change significantly over the long time spans. When the "star" passes 
through the equatorial plane above the disc's inner radius, it would be also affected by the disc mecha nically, but we do not try to m odel 



this interaction here (only gravitational effect is taken into account). This complements the approach of IVokrouhlickv & KarasI ( 119981) who 
considered Newtonian description, while taking the mechanical interaction into account. 

Let us stress, however, that we try to get a first idea about the system here, not to model any actual astrophysical situation. Namely, in 
scanning through the phase space, the ranges of parameters are not necessarily limited to "realistic" values. 



© RAS, MNRAS 000,[T]-?? 



Geodesic chaos around black holes with discs or rings 1 1 



5 POINCARE SECTIONS 

It can be expected that the originally regular, Schwarzschild geodesic dynamics will gradually grow chaotic when increasing the external- 
source (ring or disc) mass. We will begin the study of the system by plotting the Poincare sections. Actually, the (t,(f>;u^ ,u'^) phase subspace 
being completely regular due to the conservation of Ut and u^, the motion is only "non-trivial" in the {p,z;u'' ,u'^) subspace. Because of the 
four- velocity normalisation 

1 a/3 -2v c'2 , 2v , 2\-2v \/ p\2 . r z\2~\ 

— 1=5 '^UaUp = — e i+e ^+'5 |.(" ) + ) J ' 

it is confined to a three-dimensional hypersurface there. This corresponds to a system with 2D Hamiltonian. In such a case, the phase-space 
layers with regular and stochastic motion neighbour with each other (along KAM surfaces) and Poincare sections are appropriate and well 
tried method of how to study their arrangement. (In higher-dimensional systems the stochastic layers form a dense web that approaches 
arbitrarily close to any regular orbit.) 

In a stationary, axially symmetric and reflectionally symmetric field, it is natural to choose the equatorial plane {z = 0) as section on 
which the passages of orbits will be recorded (namely the radius and the radial component of four-velocity are usually registered). Orbits 
with I = Q may also be registered on the symmetry axis. Each individual Poincare section will show passages of a set of particles with given 
energy and azimuthal angular momentum (given £ and £)0 

The existence and nature of the constants of motion imply that certain subsets of geodesies will remain regular, whether the additional 
source is present or not. Specifically, £ is purely connected with it' and (. with , plus there is the four-velocity normalisation, so the 
regularity applies to the geodesies lying fully within the 3-dimensional world-sheets spanned by t, 4> and one of the meridional coordinates. 
For our systems this factually means the geodesies bound to the equatorial plane (z = 0, it^ = 0): due to the reflectional symmetry, = Q 
is their fourth integral of motion. 

We will study, on Poincare sections, the change of time-like geodesic structure in dependence (i) on relative mass and position of the 
external ring/disc and (ii) on constants of the geodesic motion (namely energy; the dependence on angul ar momentum is less intriguing). 
The geodesies are computed using the (adjusted) programme written by M. Zacek for our earlier paper l lSemeraketai]|l999l) where the 



motion of free test particles was followed in similar fields we deal with here. In all figures, a set of several hundreds of particles is launched 
from a certain region of relevant radii in the equatorial plane, with given £ and t (which determine u* and u**, respectively), with several 
different values of within the local light cone and components fixed accordingly from the four-velocity normalisation. In most plots, 
the particles are followed for about ISOOOOAf — 250000i\f of proper time; this usually corresponds to some 200 — 1000 orbital periods, i.e. 
to 400 — 2000 recorded points (since our section 2 = is the plane of symmetry, we can register passages in both directions). 

Let us note that it is no point to specify the initial conditions and evolution of individual orbits precisely, because the figures cannot be 
compared "invariantly" anyway: they correspond to different space-time backgrounds (this applies to sequences showing dependence on M 
and on rin) or at least to different regions of permitted u'^{r) (as in sequences showing dependence on £). This also implies that the plots 
naturally do not represent the same number of particles and neither the same average "density" of crossings. In particular, those with the 
permitted (r,M'^)-region opened towards the horizon tend to be more "diluted" due to the escape of particles into a black hole. 

The geodesic dynamics of perturbed systems is usually rather complicated and one has to look through many Poincare sections in order 
to catch its evolution with the parameters. Even single sections are often quite rich and require enough space in order to be represented 
properly. We hope the following groupings of plots are a reasonable compromise between providing accurate information and keeping the 
paper tolerably long. 

Figures |4l |5] and [6] concern the space-time of a Schwarzschild black hole surrounded by the inverted first Morgan-Morgan disc; they 
illustrate the dependence of geodesic dynamics on relative mass of the disc M/A/, on its inner radius rdisc and on energy of the geodesies 
£, respectively. For comparison with different discs, figures |7] and [8] show the dependence on M/M and rdisc for the Schwarzschild black 
hole surrounded by the inverted fourth Morgan-Morgan disc, and figure |9] shows the dependence on M/Af for the hole with the (n=4, m=4) 
power-law-density disc. Finally, figures [TO] [TT] and [12] provide the same kind of information for fields involving the most concentrated, 
singular external source — the Bach-Weyl ring. In all the figures, the Schwarzschild radius r is given along the horizontal axis in units of the 
black-hole mass M. The boundary of the accessible region of phase space is indicated in red. 

The presented sequences of plots can be best annotated by quotation from the introduction of the book Lichtenberg & LiebermanI ( 1992 ). 



When describing the phase-space portrait of nearly-integrable systems and KAM theory, the authors write: "The invariant surfaces break their 
topology near resonances to form island chains. Within these islands the topology is again broken to form yet other chains of still fines islands. 
On ever finer scale, one sees islands within islands. But this structure is only part of the picture, for densely interwoven within these invariant 



^ When deciding the way of choosing the initial conditions, the question always arises of what trajectories "go together". This surely depends on the problem 
solved, but one usually has two major possibilities of how to define the given set of trajectories — by fixing either their "local", or "global" parameters. When 
studying space-time motion, the "local" choice (launching the particles with given velocity in various directions, as measured by a local physical observer) 
usually feels more physical, because it more closely resembles the actual astrophysical processes and does not rely on quantities referring to conditions that 
might nowhere occur along the trajectories (e.g. like referring to "infinity"). However, the problem of regular/chaotic dynamics is tightly connected with the 
existence and values of integrals of motion, so we will adhere to these "global" quantities in assigning the initial conditions. Fixing such "global" paraineters 
also more probably confines the orbits to similar regions in phase space, which makes the comparison of their long-term evolution easier. 
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Figure 4. (The figure is on previous page!) Geodesic dynamics in the field of a Schwarzschild blaclc hole surrounded by the inverted 1st Morgan-Morgan disc 
(with radius r^jgc = 20M): dependence on relative disc mass Ai/M (its value is indicated in the plots). Passages of orbits with £ = 3.75Af, £ = 0.955 
through the equatorial plane are drawn. The boundary of the accessible region of phase space is indicated in red (this appHes to all figures). 



Figure 5. (This page.) Geodesic dynamics in the field of a Schwarzschild black hole surrounded by the inverted 1st Morgan-Morgan disc (with mass A4 = 
0.5M): dependence on disc inner radius r^n^c (its value is indicated in the plots). Passages of orbits with £ = 3.75M, £ = 0.955 through the equatorial plane 
are drawn. The figure continues on the next page. 
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Structures are thin phase space layers in which the motion is stochastic." (It should be remarked, however, that the KAM theory usually 
assumes continuous force field, which is not satisfied along orbits crossing the disc.) 

We do not present any systematic comparison of different additional sources nor a thorough account of dependences on parameters here, 
the selected figures should just illustrate overall tendencies. Although they clearly contain a number of interesting details, it would be too 
lengthy to take notice of them individually. Therefore, let us only mention rather generic observations: 

(i) Dependence on angular momentum. Larger angular momentum means larger part of (some given) energy allotted to azimuthal motion. 
This component of motion is "held" by the conserved value of £, however, so larger angular momentum favours regularity. Larger £ also 
means higher centrifugal barrier around the black hole; for a certain value of £ the barrier even prevents the particles from falling into a hole. 
This is advantageous for the present study since the particles then do not "drop out" of the phase space and can be followed for a sufficient 
number of orbital periods. (The dependence on £ is otherwise not so interesting and we do not focus on it. When exploring the dependence 
on other parameters, we mostly choose £ — 3.75M.) 

(ii) Dependence on energy. The picture begins with one regular island around an orbit with periodicity k = 1. With increasing energy, 
the primary island enlarges and shifts to higher radii, smaller islands of higher-periodic orbits arise within it and stochastic layers gradually 
develop from separatrices around them (mainly from the letter's vertices, i.e. hyperbolic points of the KAM theory, and near the boundary of 
the accessible region first). Stochastic layers tend to coalesce and form a "chaotic sea" while smaller islands of higher and higher resonances 
emerge (and often disappear again). At certain phase the picture is rather complicated, containing both the residues of regular islands of 
various sizes, narrow stochastic layers of "almost regular" trajectories (which remain very close to some regular orbit for a long time, while 
only slightly changing their action), and also large "chaotic sea" of orbits with strongly variable action (usually spread from the boundaries 
of the phase-space region accessible to the particles with given constants of motion). For large energies the higher-resonance islands continue 
to shrink towards the margin of the accessible zone, but the primary island begins to grow back and finally occupies almost all phase space 
again. This is also connected with the evolution of the accessible region (its boundaries are indicated in red in the figures). For small energies 
it typically consists of two parts, one just above the horizon and the second below (and/or around) the external source. With increasing energy, 
the two lobes get closer and finally join, permitting the particles to fall below the horizon. This effectively drains the chaotic particles (the 
most eccentric ones) from the system. 

(iii) Dependence on mass of the additional source (strength of the perturbation). The behaviour is similar to the dependence on energy: 
with increasing relative mass, higher and higher resonances gradually appear, release into a developing chaotic sea and shrink towards the 
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Figure 6. Geodesic dynamics in the field of a Sciiwarzscliild black hole surrounded by the inverted 1st Morgan-Morgan disc (with Ai = 0.5M, rdisc = 
20M): dependence on orbital energy at infinity £ (its value is indicated in the plots). Passages of orbits with £ = 3.75M through the equatorial plane are 
drawn. 

boundary of the accessible zone. The degree of chaoticity is the highest when the masses of both sources are comparable. (Even the k — 1 
orbit may become unstable and its island may break up, giving way to chaos in the central region.) When the total gravity of the external 
source prevails, however, the system rather inclines back to regularity — the central "circular" orbit again turns stable and its island gradually 
spreads out. 

(iv) Dependence on position of the additional source. The position of the ring/disc may affect the total potential even more than the 
external-source mass and this also reflects on geodesies. When the source is very close to the hole, it does not "develop" its own potential 
minimum with "its own" periodic orbits (unless it is extremely massive), so the dynamics remains almost regular (the particles only feel the 
centre heavier). Shifting the source farther from the hole, the potential valleys with periodic orbits originate. Their islands and separatrices 
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Figure 7. Geodesic dynamics in the field of a Scliwarzschild black hole surrounded by the inverted 4th Morgan-Morgan disc (with rdisc = 20M): dependence 
on relative disc mass M /M (its value is indicated in the plots). Passages of orbits with £ = 3.75M, £ = 0.968 through the equatorial plane are drawn. 

evolve in a complicated manner, making way to stochastic layers and arising again elsewhere. However, when the ring/disc is above some 
radius, the chaotic orbits begin to gather into narrow belts again, invariant curves finally emerge and form a primary and the higher islands, 
gradually spreading towards the margins of the accessible region. 

Several more general points: 

(i) More compact source perturbs the dynamics more. The Poincare-section series thus reveal slower evolution with parameters for sources 
that are more "diffuse". Of the discs considered here, the inverted 1st Morgan-Morgan disc is more compact than the 4th, the power-law disc 
(n = 4:,m — 4) being quite similar to the latter. The Bach-Weyl ring is of course the most compact perturbation and deserves a separate 
comment. The ring being singular, it always has a certain potential valley of its own where particles may evolve without passing over to 
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Figure 8. Geodesic dynamics in the field of a Schwarzschild black hole surrounded by the inverted 4th Morgan-Morgan disc (with Ai = 0.5M): dependence 
on disc inner radius r^isc (its value is indicated in the plots). Passages of orbits with £ = 3.75M, E = 0.964 through the equatorial plane are drawn. 

the black-hole sphere of influence. In this region surrounding the ring the phase-space structures behave similarly as we described above, in 
particular, the geodesic motion may become chaotic there to a large measure. At the same time, a completely regular closed region may exist 
around a potential minimum between the hole and the ring — see e.g. figure[TT] (In our figure series the accessible lobes are mostly joined, 
however.) The figures show that the hole-ring phase portrait is extremely rich and very strongly dependent on the parameters. Although the 
overall tendencies are similar to those observed with discs, there are much more and distinct higher resonances, already occurring at lower 
values of the parameters. For instance, the ring of mass M — 0.001i\f already induces a considerable perturbation of the phase pattern in 
almost the whole volume of the accessible zone (cf. figure [TOt. However, no chaotic sea is formed until values around A4 = 0.1 AI, the 
numerous regular islands and chaotic layers remain densely interwoven. 

(ii) The dynamics does not respond to the change of perturbation parameters in a "monotonous" (and by no way "linear") way. Within 
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some intervals of the parameter space just minor quantitative changes occur, whereas elsewhere the picture "quickly" alters considerably (reg- 
ularity islands suddenly disappear, sometimes arising in a totally different arrangement in a moment, indicating change in orbital periodicity. 
(A similar remark may even apply to one single trajectory in a one single background — see section[6]below). 

(iii) When increasing a given parameter, one often notices several periods of prevailing regularity and several rather stochastic ones; 
however, the dynamics typically incli nes to regularity for b oth very small and very large values of the parameters. (This behaviour has also 
been observed in other systems — cf. Strans kv et al. 1 200^, for example.) 

(iv) For our gravitational systems, the geodesic dynamics rather tends to "break up" from the boundaries of the accessible phase-space 
region, while a certain regular region mostly (though not always) survives in the interior. 

(v) The Poincare sections are generally not symmetrical with respect to u'' = 0. This specifically applies to those cases where the 
accessible region of phase space is open towards the black hole. Namely the asymmetry is due to the particles that fall into the hole (notice, in 
those cases, the excess of dots with ii' < at low radii); these are not necessarily balanced by counterparts that would be launched outwards 
from near the horizon. 

(vi) We also drew several axial sections 6 = Q for orbits with zero angular momentum, but these are not presented here. 



5.1 A comment on numerics 

We solve the geodesic-equation system (I14b~lll7b by the Runge-Kutta (actually the Hu(a) 6th-order method with variable proper- time step, 
on the cluster Tiger operated at the Astronomical Institute of our university. The numerical and computational demands strongly depend on 
specific space-time background. The case with Bach-Weyl ring and with 1st inverted Morgan-Morgan disc can be integrated quite fast and 
with high precision — constants of the motion and four- velocity normalisation are conserved there with relative error 10"^'' or less. For the 
4th inverted Morgan-Morgan disc the error is about an order larger. For the power-law-density disc the computation takes about 30 times 
longer than with the previous sources and the error is larger, about 10^^" (which is not always sufficient as seen in figures|9] contours within 
the regular regions often get blurred there at places where the phase space is more "dense"). This is because its potential is given by an 
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Figure 10. Geodesic dynamics in tlie field of a Schwarzscliild black hole sun'ounded by the Bach-Weyl ring (with radius r^-^g = 20M); dependence on 
relative ring mass M/M (its value is indicated in the plots). Passages of orbits with £ = 3.75M, £ = 0.977 through the equatorial plane are drawn. The 
figure continues on the next page. 



infinite sum of Legendre polynomials, namely l ll2b or U3b : the latter converge rather slowly and have to be cut somewhere. In the code we 
add up all terms which are greater then 10"^'^ in both the sums. Here we mean the sums themselves, without the coefficients in front. (For 

(1+-) 

the particular disc we consider, i.e. n = 4, m = 4, the front coefficient — ^ "* ^ is approximately 9.7A4 /b which is typically of the order 
0.01 ~ 1.) For comparison: at place of the worst behaviour, y^fp + z^ ~ b, the magnitude of the leading (j — 0) term of the sums is about 
0.034. Trying to estimate the tendency of farther terms, one finds that on the axis (p = 0) at |z| = b, the j-th term of the sum l ll2b reads 

(^^^) +1 ^ 2^^3' (io^^' ; for n = 4, m = 4 this falls off as ~ at large values of j. (We note in passing that 388 terms have to be 
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summed there according to our 10~ criterion; close to the disc, even more of them may be necessary. On the other hand, in regions farther 
from the \J ^ ~ h interface, several tens of terms usually suffice.) 

Finally, let us add that the numerical integration is not accurate close to the horizon (regardless the type of the external source): notice 
that some of the particles plunging into the hole got beyond the boundary of the accessible zone. 



6 TIME SERIES AND THEIR FOURIER SPECTRA 

"Chaoticity" of the system may surface in many respects. Its various characteristics tend to evolve in a more "random" way when they are 
not bound by a sufficient number of isolating integrals. One of simple tests is thus to plot the time series of selected quantities and their 
Fourier spectra. For "chaotic" systems, one expects the series to exhibit more randomness and the spectra to contain more "lower harmonics" 
(possibly even blurred clusters of frequencies). We have picked the radial velocity and vertical position, and also the particle's "latitudinal 



© RAS, MNRAS 000,[T]-?? 




5 10 15 5 10 15 20 25 



Figure 11. Geodesic dynamics in the field of a Schwarzschild black hole surrounded by the Bach-Weyl ring (with mass Ai = 0.5M): dependence on ring 
radius rri^g (its value is indicated in the plots). Passages of orbits with £ = 3.75M, £ = 0.94 through the equatorial plane are drawn. The figure continues 
on the next page. 

action", obtained by averaging the latitudinal four- velocity over an orbital period. (For these purposes, we followed the geodesies for a longer 
span, up to a coordinate time of about 10® M which typically corresponds to several thousands orbital periods.) 

The easiest time series are those of phase variables. We have plotted the vertical position z and the radial velocity plus their Fourier 
spectra for several orbits encountered in the previous (Poincare) section. (We show only the z behaviour here, but the ii'' series behave 
similarly.) In Fig.[T3] we pick two orbits around a black hole with the inverted 1st Morgan-Morgan disc and plot their spatial traces, Poincare 
sections, time series of their z position and the letter's Fourier spectra. For the regular orbit (which produces smooth curves on Poincare 
surfaces), the time series follows regular oscillations, its Fourier spectrum only contains several distinctive peaks and goes over to nearly 
constant function at low frequencies. Chaotic orbit (which produces densely dotted, "stochastic" layers on Poincare surfaces), on the other 
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hand, corresponds to complicated, inharmonic time series and fuzzy Fourier spectrum spanning over "all" frequencies including the low 
ones. The irregularity of evolution can be seen at various scales, it need not be only "local" — from time to time, the orbital parameters may 
"suddenly" change substantially and then calm down again for a longer period. This type of behaviour is very important for the evolution of 
astr ophysical systems (it i s known from the backtracking of the Solar system history, for example). 



Kovama et al.l j2007f) found, in studying the dynamics of spin particles in the Schwarzschild field, that two cases can be further distin- 



guished within irregular orbits according to the low-frequency shape of their time-series spectra. The "highly chaotic" orbits (filling large 
area — the "chaotic sea" — in Poincare sections) show "white noise" there, i.e. the noise with relatively low and approximately constant 
mean value. On the other hand, the "less chaotic" orbits (filling just narrow chaotic layers, namely confined to the vicinity of regular islands 
— this case is being called "sticky motion") show the 1// dependence at low frequencies (the spectrum begins at higher amplitudes and goes 
down with frequency). We have checked whether similar behaviours can be found in our systems, too. 

In figure [14] we compare several rather low-energy {£ = 0.953) orbits from the "moderately chaotic" phase space of a black hole 
surrounded by the inverted 1st Morgan-Morgan disc (M — 0.5Af, raise ~ 20M). First we took "most regular" orbits that live deep in the 
primary island, almost unaffected by perturbation (pink, light blue and violet). Their spectrum is very smooth, containing just a few peaks at 
nearly the same frequencies. The second triple (orange, blue and green) comes from the regular island of periodicity k — 3. The spectra are 
again smooth with distinct peaks at similar frequencies, now containing more of them (these orbits only exist due to the perturbation and are 
more complicated then those belonging to the primary island). "Even less regular" are the trajectories living on the outskirts of the primary 
island, encircling also the higher-periodic (e.g. the fc = 3) orbits and neighbouring with the chaotic layer (These trajectories break up into 
the chaotic sea when energy is increased.) Their spectra already contain much more harmonics (near the separatrices, the resonances are very 
close to each other), but the low-frequency part is still smooth and constant. Finally, we took three orbits from the stochastic layer which, 
however, keep close to the k = 2 islands and so might be only "weakly chaotic" (black, red and green on the second Poincare section on 
the right; mainly the black trajectory is seen to linger in the vicinity of the regular islands for a long time). The spectra of all these orbits are 
clearly chaotic and show the 1// depende nce at low frequencie s. 

The figures confirm that the results of lKovama et al J 1 120071) apply to our systems as well. We wanted to see in further detail the difference 



between the "almost regular" black trajectory, sticking to the k — 2 island for a considerable amount of time, and the truly regular trajectory 
belonging to that island. In figure[T5] one can see clear difference between time series of z obtained for these seemingly very similar orbits. 
However, in proper times between r — 6.7 ■ 10^ Af and r — 9.4 • lO^M the chaotic orbit really behaves very much like a regular one. This is 
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Figure 12. Geodesic dynamics in tlie field of a Schwarzscliild black hole sun'ounded by the Bach-Weyl ring (with M = 0.5M, r^ing = 20M): dependence 
on orbital energy at infinity S (its value is indicated in the plots). Passages of orbits with £ = 3.75M through the equatorial plane are drawn. The figure 
continues on the next page. 



also confirmed by spatial picture of the orbits in the figure: in the above stage, the chaotic orbit is almost indistinguishable from the regular 
one, although it clearly lo oks chaotic in other times. 



Kovama et al.l ( I2007D suggested that the degree of orbital chaoticity depends on strength of the perturbation (magnitude of the particles' 



spin in their case) rather than on energy. In order to check this on our system, we compared a few more orbits computed in the field with 
much heavier disc (A4 — 1.3A/) and with constants of motion for which the phase space is strongly chaotic. In figure[T6] eight orbits filling 
the chaotic sea are compared with four orbits spending most of the time in the vicinity of regular islands. Power spectra of the ^-position time 
series are different and indicate that the first set belongs to the white-noise type, whereas the second set belongs to the 1 // type. Hence, even 
in a more strongly perturbed and thus rather chaotic system there exist trajectories that "mimic" the behaviour of regular ones for a certain 
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amount of time (the action changes only shghtly there), and this reflects in the 1 // shape of their low-frequency power spectra, even if they 
may inhabit larger portions of phase space at other times. The power- spectrum behaviour is thus a property of the phase-space layer to which 
the given orbit belongs, rather than of the space-time — or the dynamical system — as such. 

Let us add that we obtained similar results for the other discs (other than the inverted 1st Morgan-Morgan), too. 



7 "LATITUDINAL ACTION" INTEGRAL 

There are various other variables whose time series should reflect the nature of the system dynamics. For example, one could register the time 
intervals between particle's returns to a given r or a giv en 9, the times between reach ing turning points, etc. In a paper on charged-particle 
motion in an EM field surrounding a rotating black hole. lTakahashi & Kovamal ( EoO^ considered the "latitudinal action" 



Jb = \ <j> ^Jgeeu^u^ dr , (27) 

given by averaging the test-particle latitudinal four-velocity over one orbital cycle (of length /). They compared its time evolution with 
Poincare maps and evolution of Q and really observed clear correlation. We were curious whether a similar observation can be made in our 
disc/ring-perturbed systems, too. 

For regular orbits, we obtained regular oscillations of the action around certain constant mean value. The amplitude of these oscillations 
is very small for orbits close to the primary resonance ("circular" periodic orbit) and larger for orbits belonging to higher-periodicity islands. 
For chaotic orbits, the action integral does not evolve in a harmonic way and is not tied to any obvious mean value. In figure[T7] the behaviour 
of action is drawn for the four "weakly chaotic" orbits from figure[T6](spectra of their z-position are given at bottom right of that figure). All 
these trajectories are "stuck on" regular islands in certain stages of their evolution — and it is seen now that in those stages their action really 
behaves almost like along regular orbits (almost regular moderate oscillations around almost constant mean value). 

Let us conclude with a trajectory that also adheres to a regular island for a certain time, but its power spectrum still does not show the 
1// shape at low frequencies. In figure[76l it is the dark blue geodesic occurring in the vicinity of the fc = 3 islands. In figure [TS] we plot the 
time series of its latitudinal action and spectrum of its z-position evolution. Clearly the action behaves quite regularly for a considerable time 
(mainly in periods from to 150A/ and from 400A'f to IBOOAf), yet the z-motion spectrum is still of clear "white-noise" type, not showing 
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Figure 13. Compari.son of a selected regular (left) and chaotic (right) orbit, both with con.stants £ = 0.955, € = 3.75M and obtained in the field of a 
Schwarzschild black hole surrounded by the inverted 1st Morgan-Morgan disc with M = 0.5M, r^^^^ = 20M. From top to bottom, the rows contain spatial 
tracks of the orbits, their Poincare sections z = (r,u^), time series of the z position and corresponding Fourier spectra. Coordinates are in the units of M, 
frequencies in the units of 1/M. 
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Figure 14. Comparison of selected orbits with constants S = 0.953, t = 3.75A/, obtained in the field of a Schwarzschild black hole surrounded by the 
inverted 1st Morgan-Morgan disc with M = 0.5M, r^isc = 20M. At top left, 3 types of regular orbits are chosen in respective Poincare section. The spectra 
of their z-position time series are given at top right (three primary-island orbits), middle left (three orbits from islands of periodicity fc = 3) and middle right 
(two orbits skirting the whole island and neighbouring with the chaotic layer). In the bottom row, three "weakly chaotic" orbits are chosen in the same Poincare 
section and their spectra shown on the right. The individual orbits are given in different colours. Radius is in the units of M, frequency in the units of 1/M. 



any signs of the 1// dependence. Perhaps the latter would only appear if the spectrum were calculated solely for the "almost regular" phase 
of the orbit. 
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Figure 15. Further comparison of the orbit sticking to the fc = 2 islands in Fig. ll4l (the black one in the bottom row) and a selected regular orbit belonging to 
these islands. On the left, the spatial track (top), the Poincare section (middle) and the time series of z (bottom) are drawn for the regular orbit. On the right, 
the spatial tracks for the weakly chaotic orbit are drawn, at "almost regular" times 7 ■ lO^M < r < 9 ■ lO^M (top) and at more chaotic times r < lOOOOOM 
(middle); the time series of z is given at the bottom. 
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Figure 16. Two types of chaotic orbits (with constants S = 0.956, £ = 4M) found in the field of a black hole surrounded by the inverted 1st Morgan-Morgan 
disc with M = 1.3M, r^^^^ = 20M. Eight orbits fining the "chaotic sea" on Poincare section (top) have their z-position power spectra at bottom left, while 
four orbits dwelling mostly in the vicinity of regular islands in the centre have their spectra at bottom right. 



8 CONCLUDING REMARKS 

We have tried to get a basic idea about how the dynamics of free massive test particles in a Schwarzschild black-hole field is perturbed by 
gravity of a static, axially symmetric and reflectionally symmetric ring or disc. Poincare sections have proved that it becomes chaotic and that 
the degree of stochasticity grows with compactness and relative mass of the external source and with energy of the particles. However, for 
large values of these parameters, the dynamics rather returns towards regular regime. The shape of stochastic regions depend very much on 
radius of the additional source. This is clearly the consequence of how the parameters affect the shape of effective potential, in particular of 
whether and where are regions that allow "permanent" orbital motion. Our system shows typical features of a nearly-integrable Hamiltonian 
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Figure 18. Time series of the latitudinal action (left) and z-position power spectrum (right) for the "dark blue" orbit from fig. [16] The orbit spends quite a long 
time very near the regular islands of periodicity k = 3 (they are not far from the boundary of the accessible region), but its power spectrum still does not show 
any signs of the 1/frequency behaviour at low frequencies. 
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system, the behaviour revealed on Poincare sections is a good illustration of the KAM theory. (Strictly speaking, the latter should only be 
applied to orbits which do not cross the disc, however.) 

The appearance of selected orbits in Poincare sections were compared with their spatial picture, with the time series representing 
evolution of their vertical (z) position, with the latter's Fourier spectra and with the evolution of the "latitudinal action" assigned to their 
individual loops (fixed by ^-periodicity). Chaotic orbits are clearly distinguishable from regular in all these respects. Within the chaotic 
trajectories, one can further recognise those which adhere to regular orbits, lingering within narrow bands in Poincare diagrams for most of 
the time, and those which fill a large area ("chaotic sea") there in a rather uniform way. The degree of chaoticity is reflected in time series 
and their Fourier spectra: the "weakly chaotic" orbits produce nearly regular and moderate oscillations around nearly constant mean value, 
which typically gives rise to "1/frequency" shape of the spectrum at low frequencies, whereas "strongly chaotic" orbits produce strongly 
irregular times series without any distinct mean value, which corresponds to "white-noise" shape at low end of the spectrum (relatively flat 
curve at rather low values). The most interesting are chaotic orbits in strongly perturbed systems which undergo both the "sticky motion" 
near regular islands and "unrestrained" phases of drifting around the chaotic sea. The power spectra of their evolution usually "take notice" 
of the almost regular phase and fall under the 1// type, but we also came across a trajectory whose spectrum showed a clear white noise, 
although it behaved quite orderly for a considerable time. 

The study of any specific dynamical system may develop in several obvious ways. First, there should be a running check of gener- 
icity/stability of the dynamical regime found. This involves going through the initial conditions in an exhaustive manner, changing the 
parameters of the system, trying different numerical schemes. . . Second, one may elaborate the model of the explored system in order to 
make it as realistic as possible, try to estimate and examine the consequences and mutual couplings of different "perturbations"/neglects 
as well as their astrophysical importance. Finally, it is very appropriate to study the dynamics of the specific system by more methods and 
observe the relations between different aspects of its chaoticity. This is our next plan. 

Concerning the last point(s), let us note that there especially arises the question about the role of features which are specific for the 
theory that underlies the phase space. In general relativity, suc h a special fe ature is the space-time curvature. It is clear that chaos is not 
solely caused by space-time curvature. (See the quotation from lVieira & Lete lier (1996b) in Sect. 12.21 As a matter of fact, some systems 
are chaotic both in the Newton's theory, in special relativity and in general relativity.) However, the example of the Einstein's theory of 
gravitation itself suggests that it might be possible to interpret the dynamics of some systems as geodesic flows in suitably defined mani- 
folds. This would certainly depend on the type of interaction, but also on some integral propertie s of the spec i fic traiectory (like constants 
of the motion). This direction has been notably pursued by Szydlowski and coworkers, see e.g. Szvdtowskil 1 1998L 1999h and references 
therein, while dynamics of specific systems were interpreted "geometrically" in some of the papers already cited in Sect. 12.21 j Yurtsevei 
1995 : Sota et al. 19 96) (cf. section 2.3 of , Cor nish & Gibbons ( 1997 ) ). More recen t contributions t o this line of research include mainly 



Casetti et all j200(]|):lRamasubramanian & Srirairj ( i200ll) Jciementi & Pettiiiil j2002h : [Kawabel j20oi) : lHorwitz et alj J2007l) : lcafard Jiooj) : 



Zion & Horwitzl ( 120081) : ICerruti-Sola et alj hoOSh (further references are g ven therein). 
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